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ABSTRACT 

The narrow GD-1 stream of stars, spanning 60 deg on the sky at a distance of ~ 10 kpc from the Sun 
and 15 kpc from the Galactic center, is presumed to be debris from a tidally disrupted star cluster 
that traces out a test-particle orbit in the Milky Way halo. We combine SDSS photometry, USNO-B 
astrometry, and SDSS and Calar Alto spectroscopy to construct a complete, empirical 6-dimensional 
phase-space map of the stream. We find that an eccentric orbit in a flattened isothermal potential 
describes this phase-space map well. Even after marginalizing over the stream orbital parameters and 
the distance from the Sun to the Galactic center, the orbital fit to GD-1 places strong constraints on the 
circular velocity at the Sun's radius V c — 224 ± 13 km/s and total potential flattening g$ = 0.87+QQ4. 
When we drop any informative priors on V c the GD-1 constraint becomes V c — 221 ± 18 km/s. Our 
6-D map of GD-1 therefore yields the best current constraint on V c and the only strong constraint on 
g$ at Galactocentric radii near R^ 15 kpc. Much, if not all, of the total potential flattening may be 
attributed to the mass in the stellar disk, so the GD-1 constraints on the flattening of the halo itself 
are weak: q^>^haio > 0.89 at 90% confidence. The greatest uncertainty in the 6-D map and the orbital 
analysis stems from the photometric distances, which will be obviated by Gaia. 

Subject headings: Galaxy: fundamental parameters - Galaxy: halo - Galaxy: kinematics and dynamics 
- surveys - stars: kinematics - stellar dynamics 



1. INTRODUCTION 

The Sloan Digital Sky Survey (SDSS) is an imag- 
ing and spectroscopy survey which mapped quar- 
ter of the sky near the North Galactic Cap. The 
data have proven extremely useful for the under- 
standing of the Milky Way halo. In addition to 
a large list of MW s a tellites (|Belokurov et al.l l2007ct 
Koposov etafl I2007ah llrwin et al.1 120071 : iWalsh et al.1 
2007ft several extended stellar sub-structures in the 
MW halo have been found in the SDSS data, such 
as the tidal tail of the Palomar 5 globular cluste r 
(jOdenkirchen et aJJ [200lt iGrillmair fc Dionatosll2006aD . 
the Monoceros ring (|Newberg et al.l 12002ft . two north- 
ern tidal arms of the disrupting Sagittarius galaxy 

Orphan" stream 
the Aquila 



Belokuroy et al.1 12006IL the so cal led 

Grillmairl 120061 IBelokurov et al.1 I2007bl) . 



overdensitv (jBelokurov et al.l 12007a ) and the very long 
thin s tellar strea m called GD-1 (Grillmair fc Dionatoa 
I2006bf ). Recently IGrillmairl (|2009f ) claimed the discov- 
ery of another 4 stellar streams. Streams are presumed 
to be remnants of tidally disrupted satellite galaxies and 
clusters. They provide important insights into the his- 
tory of accretion events and the physics of Galaxy forma- 
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tion. The tidal debris from disrupted satellites (clusters) 
spreads out in orbital phase on a path that is close to the 
orbit of the progenitor. Streams tracing out orbits there- 
fore provide opportunities to constrain the Milky Way's 
gravitational potential. 

After in itial searches for tida l tails of globular clus- 
ters (e.g. IGrillmair et al.l [1995) it was the extended 
Sagittarius tidal tail that first made deriving such con- 
straints practical (see e.g.llbata et al.H2001t iHelmil [20041 : 
Uohnston et al.l l2005t lLaw et al.l l2005h . However, the 
tidal tail of the Sagittarius galaxy is quite wide and con- 
tains a considerable mixture of different stellar orbits, 
making it complex to model. For constraining the gravi- 
tational potential, a stellar stream that is very thin but of 
large angular extent, is ideal, because it permits precise 
orbital models. 

The first studies of globular cluster tidal debris only 
revealed short (< 1°) signs of tails, but in recent years 
with the advent of large photometric surveys such as 
SDSS and 2MASS and significant advances in the tech- 
niques used to find streams, significan t progress has been 
made. The matched filte r technique (Odcnkir chen et all 
l200lHRockosi et al.ll2002T > has revealed the beautiful tidal 
stream of Palomar 5. Detail ed analysis of the Pal 5 
stream, inc luding kinematics (|Odenkirchen et all I2001L 
I2003L 12009ft . have shown the promise of this approach, 
but also revealed that data over more than 10° on the 
sky are needed to pla c e good constraints on t he pot en- 
tial. IGrillmairl (1200 6). Gri llmair fc Dionatosl (|2006al |bfj . 
IGrillmair fc Johnson! (|2006l ) were successful in the detec- 
tion of very long stellar streams using this technique, 
including the 63° long stellar stream GD-1 . Besides 
the stream length and the approximate distance, most of 
the properties of GD-1 were unknown. Since the stream 
is long but relatively thin, with no apparent progenitor 
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remnant, it was suggested that it arose from a globular 
cluster. In this paper we make an attempt to determine 
all possible properties of the GD-1 stream including dis- 
tance, position on the sky, proper motion, and radial 
velocity and try to constrain the Milky Way potential 
using that inform ation. This work goe s in parallel with 
the work done bv lWillett et al.l (|2009t ). but we are able 
to get a full 6-D phase space map of the stream and are 
able to use that map to provide significant constraints 
on the MW potential. See also lEvre fe Binnevl (|2009l ) 
for theoretical discussion of using thin streams in order 
to constrain the MW potential. 

In performing this study we have obtained the first 6-D 
phase-space map for a kinematically cold stellar stream 
in the Milky Way. We view our p resent analysis in same 
sense as a pilot study for the Gaia (jPerrvman et al.ll2001|) 
age, when this ESA space mission will deliver dramati- 
cally better data on streams such as GD-1. 

This paper is organized as follows: In Section[5]we dis- 
cuss the analysis of the SDSS photometry, which entails 
mapping the GD-1 stream in 3-D as well as determin- 
ing its stellar population properties. In Section [3] we 
present the kinematics, with proper motions from SDSS- 
USNOB1.0 and line-of-sight velocities from SDSS and 
Calar Alto. In Section U we combine this information in 
an iterative step that involves improved stream member- 
ship probabilities, which in turn affects the estimates of 
proper motions and distances. This procedure results in 
the most comprehensive 6-D data set for a stellar stream 
m our Milky Way. In Section [5] we model the stream 
data by a simple orbit in a simple parametrized grav- 
itational potential. We measure the potential circular 
velocity and find that the overall Milky Way potential 
at the GD-1 stream position is somewhat flattened, but 
that much of that flattening can be attributed to the 
disk. 

2. STELLAR POPULATION OF THE STREAM 

The probability that a star is a member of the 
GD-1 stream depends on its 6-D position and its metal- 
licity. In the space of photometric observables, this 
means that it depends on (a, (5), magnitude and color. 
In practice, the determination of the stream's angular 
position, distance and metallicity (presuming it is 'old') 
is an iterative process which we detail here. 

IGrillmair fe Dionatosl (2006b) made the initial map 
of the stream using a matched color-magnitude filter 
based on the CMD of M13 observed in the same fil- 
ters. Not presuming a particular metallicity (e.g. that 
of M13), we start our analysis with a simple color- 
magnitude box selection for stars (0.15 < g—r < 0.41 and 
18.1 < r < 19.85). The resulting distribution is shown 
in Fig. [TJ That particular color-magnitude box was se- 
lected as appropriate to find metal-poor main sequence 
(MS) stars at a distance of ~ 10 kpc, and indeed the 
stream is marginally discernible in the Figure. With just 
a color-magnitude box, however, the detection fidelity of 
that stream is notice a bly low er than that achieved by 
IGrillmair fc Dionatosl (j2006blU their Fig. 1). The distri- 
bution of stars on Figure [TJ is plotted in a rotated spher- 
ical coordinate system (fa, fa), approximately aligned 
with the stream, where fa is longitude and fa is the 
latitude. The north pole of that coordinate system is lo- 
cated at a p =34°.5987, S p = 29°. 7331, the zero-point for 
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Fig. 1.— The number density of SDSS DR7 stars with 0.15 < 
g — r < 0.41 and 18.1 < r < 19.85, shown in the rotated spher- 
ical coordinate system that is approximately aligned with the 
GD-1 stream. The map was convolved with a circular Gaussian 
with a = 0.2°. The gray arrows point to the stream, which is 
barely visible in this representation, extending horizontally near 
fa = 0° , between fa = -60° and 0° . 

fa is located at a = 200°, and we will use this coordinate 
system for convenience throughout the paper to describe 
stream positions (the transformation matrix from (a, S) 
to (fa, (j>2) is given in the appendix). 

If we integrate the low-contrast 2D map in Fig.QJalong 
the fa axis, creating a one-dimensional profile of the 
stream, the presence of the stream becomes very clear. 
Figure [2] shows this profile for stars with 0.15 < g — r < 
0.41, 18.1 < r < 19.85 and with -60° < fa < -10°. 
In that Figure we also overplot the Gaussian fit to this 
profile with ~ 600 stars and Gaussian width (a^ ) of ~ 
12'. This number of stars corresponds to a total stellar 
mass of A/* w2x 10 4 M Q , if we assume a distance of ~ 
10 kpc (see below), and a Chabrier IMF (|Chabrierll20"0lh 
with an old, metal-poor stellar population. Given that 
number of stars, we expect to see around 3000 stream 
stars in SDSS with r <22. The mean surface brightness 
of the stream is around 29 mag/sq.arcsec. 

We expand this approach to the determination of the 
CMD of the stream. The data and the fit shown in Fig- 
ure [2] was obtained for a fairly wide color- magnitude se- 
lection box. But we can construct such a profile for any 
other (e.g.) color magnitude box and that profile can 
then be fitted by 

N obs (fa\CMD) = N BG (CMD) + 
N stream (CMD) x -=L exp ( -0.5 f ^'^A \) 

V27TO-02 \ V J J 

, where CMD refers to a given color-magnitude bin, and 
where we assume that both center (fa,o) and width (cr^) 
of the stream are fixed at and 12 arcminutes. A fit of 
the Eq.Q]model to the observed data N bs(fa\CMD) can 
be performed in x 2 sense. As a result N strea m{CMD), 
the number of stream stars (and its error), can be deter- 
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Fig. 2. — One-dimensional stellar density profile across the 
stream using the stars with 0.15 < g — r < 0.41 18.1 < r < 19.85 
across the 02 = 0° axis, integrated along the stream in the interval 
-60° < 0i < -10°. The Gaussian fit with ~ 600 stars and cr=12' 
is shown in red. 

mined for each given color-magnitude bin, resulting in a 
Hess diagram for different pairs of SDSS filters (it — g, 
g — r,r — i, i — z). Figure [3] shows the resulting Hess 
diagram of the stream derived in several bands. These 
clearly show a main sequence(MS). The location of the 
MS turn-off cannot be clearly identified, although there 
may be a hint at q = 18.5, u — q = 1. On Figure [3] 
we also overplot the lMarigo et al.l (|2008l ) isochrones with 
age= 9 Gyrs, log(Z/ Z w ) = —1.4 at 8.5 kpc, which seem 
to match quite well. Ilvezic et al.l (|2008f) recently showed 
that the location of MS stars in the (it — g) — (g — r) color- 
color plane is a good metallicity diagnostic. Therefore, 
we construct the (u — g) — (g — r) color-color diagram of 
the stream stellar population shown in Fig. |4j which ex- 
hibits a distinct concentration of stars at (0.8,0.3). This 
argues for a population of single or a dominant metal- 
licity and we can conv ert this color locat ion to a metal- 
licity using Eq. 4 from Ilvezic et al.l (j2008D : [Fe/R] phot = 
— 1.9 ± 0.1. This provides a metallicity estimate that is 
directly linked to SDSS spectral metallicity estimate. 

To derive the metallicity, age, and distance of stream 
stars in a systematic way, we fit the color-magnitude 
diagrams using a grid of i sochrones populated realisti- 
cally according to the IMF (|DolphiiJl200alde Jong et al.l 
120081) . We focus on fits to the color-magnitude diagrams 
in it, g and g, r filters, since that the u — q color of the 
MS t urn-off is a good metallicity indicator (jlvezic et al.1 
2008) . We create the synthetic Hess diagrams for a 
grid of model stel lar populations (|Girardi et al.l 120001 : 
IMarigo et all I2008D 6 with different ages (3 - 12 Gyr), 
metallicities (Z = 0.0001-0.025), distances (6-14 kpc), 

6 To retrieve the isochrones we used the web interface pro- 
vided by Leo Girardi at the Astronomical Observatory of Padua 
http : //stev. oapd. inaf . it/cgi-bin/cmd_2 . 1 



and a Chabrier IMF (|Chabrierl 12001). We then explore 
that grid by computing log-likelihood of the distribution 
of stars in color-magnitude space. Figure [5] shows the 
2D profile likelihoods contours of the age vs metallicity, 
age vs distance and distance vs metallicity planes. The 
filled circle indicates the best fit model: age= 9 Gyrs, 
log(Z/Z Q ) = —1.4 and distancc= 8 kpc. Clearly the 
age is the least well constrained parameter; the distance 
seems to be relatively well constrained, but has a co- 
variance with [Fe/H]. We will revisit this issue later, as 
the analysis of Fig. Q] implies a lower metallicity. Fig. [3] 
shows that the isochrones are reproducing the observed 
Hess diagrams well, and hence further in the paper we 
will use t= 9 Gyrs, log(Z/Z Q ) = —1.4 as the baseline 
model for the stream's stellar population. It should 
be noted that the distance measurement from Fig. [5] 
represents the averaged distance along the stream from 
—50° < (j>i < —20°. In section 0] we will present esti- 
mates of the distance to different parts of the stream. 

It is noticeable that the metallicity derived from the 
CMD fitting is higher th an from the estima te based on 
empirical calibration of Ilvezic et al.l (|2008l) (see above) 
and higher than the measurement b ased on the SDSS 
spectra given bv lWillett et al.l J2009). This discrepancy 
is understandable given the known inaccurac ies of the 
isoch rones in the SDSS photometr ic system (lAn et al.l 
120081) . In particular the Figure 19 of I An et al.l (|2008l ) pa- 
per clearly shows a mismatch between fiducial isochrone 
derived for the M92 globular cluster (which is used 
elsewhere as a good approximation of old metal-poor 
stars in the halo) and the theoretical isochrones. For 
main sequence stars below the turn-off (which are of 
the most interest here), the mismatch between a fidu- 
cial isochrone of the M53 glo bular clu s ter (w hich has 

metallicity [Fe/H] 2) from I An etail (|2008t l and the 

isochrone, which we are using can be approximated by 
a distance shift of ~ 10%. Therefore in our analysis 
we reduce all the distances derived on the basis of the 
CMD fit by 10%. Due the these described inaccuracies 
with the isochrones in SDSS filters, some remaining sys- 
tematic error i n dist ances may still exist, although the 
study of lEvrel (|2009D seems to indicate it is small. Until 
isochrones in SDSS filters are fixed, the usage of fiducial 
isochrones may give more consistent results, but in the 
current paper we decided to proceed with the theoretical 
isochrones and 10% distance correction. 

Splitting the CMD data into 0i bins shows that there 
is a distance gradient along the stream: Figure [5] shows 
two Hess diagrams obtained for two different pieces of 
the stream, on the left —40° < 0i < —20°, and on 
the right -10° < 0i < 10°. The baseline isochrone is 
shifted to distances of 8.5 kpc(left) and 11 kpc(right) re- 
spectively. Clearly, the part of the stream at —10° < 
0i < 10° is further from the Sun, (as already n oted by 
IGrillmair & Dionatosll2006b1: IWillett et al.ll2009l ). 

The determination of the CMD properties for the 
stream allows us to select the possible stream mem- 
ber stars with much less background contamination, 
compared to a simple color magnitude box. Figure [7J 
shows the map of the stream (in 0i , 02 coordinates) af- 
ter applying a matched filter based on the CMD of the 
stream. For Fig. [7J each star in the SDSS dataset was 
weighted by the ratio of the stream membership probabil- 
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Fig. 3. — Color magnitude (or Hess) diagrams of the stream derived by statistical background subtraction using the Eq.[l]fit, in different 
filters (u — g vs g, g — r vs r, r — i vs i and i — z vs z (from left to right). The grayscale shows the number of stars per rectangular bins. All 
the magnitudes are extinction corrected. Overplotted are theoretical isochrones for age = 9 Gyrs, log(Z/Zg) = —1.4, distance = 8.5 kpc. 
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Fig. 4. — (u — g) — (g — r) color-color diagram of the stream, 
which constitutes a photometric metallicity estimator (following 
Ivczic ct al. 2008, Eq. 4), shown after statistical background sub- 
traction as for Fig. \3\ All the magnitudes were extinction cor- 
rected. The grayscale shows the number of stars per bin, with 
a distinct concentration of stars at (0.8,0.35), that implies a well 
defined metallicity [Fe/H] = -1.9 ± 0.1. 

ity and the background probability f f™[„^fr^,"r^ 
, where P s tream{u — g, g — r,r — i,r) have been com- 
puted based on the stellar population fit (Fig. [5]), and 
Pbg{ u ~ 9i .g — — i , r) was computed e mpirically from 
the re gions adjacent to the stream (see e.g. lRockosi et al.l 
120021 for the application of similar method) . The result- 
ing image after the CMD weighting shows the stream 
with obviously greater contrast than Figure [TJ It is no- 
ticeable that the density of stars in the stream does not 
seem do be constant and even does not change monoton- 
ically along the stream, instead the stream seems to have 
a few clumps and holes. The nature of the substructure 
in the stream is unclear. It may be related to either the 
history of the disruption process (|Kupper et al.ll2010h . or 
e.g. interaction with dark matter subhalos around MW 
(|Carlberdl20M) . 

3. KINEMATICS 

In this section we describe how we derived estimates of 
the 3-D kinematics of the stream by looking at the proper 



motions and radial velocities of the probable stream 
members. 

3.1. Proper motion 

Despite the distance of ~ 10 kpc to the GD-1 we 
demonstrate in this section that it is possible to de- 
rive good constraints on its proper motion. The anal- 
ysis is based on proper motions der ived from combin- 
ing USNO-B1.0 (IMonet et al.l 120031) with SDSS data 
(IMunn et al.l 12004 120081) which we take from SDSS 
DR7 (|Abazajian et al.ll200l . 

At a distance of ~ 10 kpc a stream that is moving per- 
pendicular to the line sight with the velocity of 200 km/s 
(roughly a characteristic velocity in the halo) should have 
a proper motion of the order of ~ 4 mas/yr. Hence the 
expected proper motion is comparable to the precision of 
individual proper motion measurements of 3 — 4 mas/yr 
(M mm ct al. 2004). As stream member stars that are ad- 
jacent on the sky should have (nearly) identical proper 
motions, we can, however, determine statistical proper 
motions for ensembles of stars. 

We start by deriving the /i-distribution of likely mem- 
ber stars, by extending the analysis of Section [5] and us- 
ing both the angular position on the sky (specifically fa) 
and the location in CMDs for each star to evaluate its 
stream membership probability. Specifically, we grid the 
proper motion space into 'pixels', then we select the stars 
within each proper motion 'pixel' and with high mem- 
bership probabilities in the CMD (u — g,g — r,r — 
space Pstream/ Pbg > 0.1. For that sample of stars 
within each proper motion 'pixel' (/2, S/i), we determine 
N s tream(i2\4>2iC'MD) (number of stream stars) by inte- 
grating along the stream direction (f>i and fitting the re- 
sulting 4>2 distribution with the Eq. [T] The resulting 
distribution N Q i, s (fl) is shown in Fig. [8j The grayscale 
in the left panel of the Figure shows the proper motion 
distribution of fi a ,fJ-s of probable stream member stars 
(Ngtream) integrated along the stream, while the con- 
tours show the proper motion distribution of the back- 
ground stars selected with the same color-magnitude cri- 
teria (corresponding to Nbg from Eq.[T]). It is clear that 
the stream stars are on average moving differently than 
the bulk of background stars. However, the observed 
proper motions contain the reflex motion of Sun's motion 
in the Galaxy. We can account for this and then con- 
vert fi to (/i^ 1 , /^0 2 ), the proper motion along the stream 
/i^ = fiaiong and proper motion across the stream 
in the Galactic rest-frame (where (f>i , 4>2 are the stream 
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Fig. 6. — Distance variation along the stream. The CMD dia- 
grams of the stream for two different parts of the stream, —40° < 
4>l < -20° (left), -10° < 0i < 10° (right). The isochrones for 
the the best fit model log(Z/Zg) = —1.4, age= 9 Gyr were shifted 
to the distance of 8.5 kpc on the left panel and to 11 kpc on the 
right panel. Some distance variation apparent, with the stellar 
population shown on the right located at greater distances. 



Fig. 5. — Isochrone fits to the color magnitude diagrams (Fig. [6} in u, g, r, i, z bands for the stream integrated over —60° < <f>i < —10°. 
The contours show the formal 60%, 90%, 99%, 99.9% confidence regions as a function of distance, age and metallicity respectively. Filled 
circles show the location of the best goodness of fit point. 

approximately know the distance to the stream, this cor- 
rection Ufa 2)C = Ha, s — Href lex can be done. We will dis- 
cuss the consequences of the uncertainties in the Sun's 
motion and the stream differences in Section [4] 

The right panel of Fig. [8J shows the distribution of 
Ufa , [ifa of the stream stars. The contours show the 
corresponding distribution of background stars with the 
similar color-magnitudes to the stream stars. Reassur- 
ingly we see that stream stars are moving approximately 
along <f>i, i.e. along the stream orbit, an important 
plausibility check for the correctness of the proper mo- 
tion measurement. In contrast, the proper motion dis- 
tribution of the background stars after subtracting the 
proper motion due to sun's movement is centered around 
{Hfa,H>fa) — (0,0), which appears reasonable since with 
our color-magnitude selection we are selecting primarily 
the halo sta rs at distances ~ 10 kpc. Those show little 
net rotation (|Carollo et all 120071 : Kue et al.l I2008T ). The 
estimate (/J, fa) ~ —8 mas/yr also immediately implies 
that the stream is going retrograde with respect to the 
Milky Way's disk rotation. 

In Figure [9] we illustrate the proper motion variation 
along the stream. These plots, showing only Hfa{<pi) 
were derived the same way as Fig. [8j except that we did 
not integrate in </>i along the entire stream but only in 
4>i intervals. The right panel of Fig. [§] shows the dis- 
tribution of the proper motions along the stream of the 
background stars. The left panel of the Figure shows 
the distribution of proper motions of likely stream mem- 
ber stars as a function of angle along the stream (the 
proper motion due to the Sun's motion was subtracted). 
The left panel reveals a slight, but significant gradient 
in (/J.fa), of the order of 3 to 5 mas/yr. Note that the 
decrease of the proper motions towards (f>i = coincides 
with the distance increase to the stream (see Fig. [6]). 

Having determined the stream proper motions, we 
can further improve the CMD-filtered map of the 
GD-1 stream (Fig. [7]) by requiring that the proper mo- 
tions of the stars are consistent with the proper motion 
of the stream. That is shown in Fig. [TO] and discussed 
in Section [321 




Fig. 7. — Map of stream stars in the rotated coordinate system 
after applying the matched CMD filter from Section fS] The figure 
shows the histogram of stars in 4>\,<j>2, where each star has been 



weighted by the CMD weight —p B 
variation of distance was assumed 



^(u — g ,g — r,r — i,r) 



-i,r) 



The linear 



coordinates introduced in Section [5]). The proper mo- 
tion component arising from the Sun's movement in the 
Galaxy can be easily computed for each star. 



Href lex 



where V is a 3-D velocity of the sun (~ 220 km/s) and 
r is the vector from the sun towards each star. As we 



3.2. Radial velocities 



6 



20 



10 



=1 



•10 



-20 




20 -' 



10 



E I 



-10 - 



-20 




-20 



-10 



Ma. 





[mas/yr] 



10 



20 



-20 



-10 





[mas/yr] 



10 



20 



Fig. 8. — Proper motion of the stream. The left panel shows the proper motion in right ascension and declination fi a ,fig (as observed, 
e.g. no correction for the Solar motion in the Galaxy was made). The right panel shows the proper motion in the rotated coordinate 
system (4>i,<f>2) (<f>i is oriented along the stream) and after the subtraction of the proper motion due to the Sun's motion in the Galaxy 
(assuming for now V c = 220 km/s). The grayscale in each of the panels shows the number of stars per proper motion bin. Contours 
corresponding to 30,60,90 stars per bin show the the proper motion distribution for the field stars, with similar colors and magnitudes 
to the stream stars. The proper motions of the stream stars are clearly distinguishable from the proper motions of the background stars. 
The right panel confirms the fact that the stream stars are moving approximately along its orbit (/10 2 ~ mas/yr), while the mean proper 
motions of background stars after subtracting Sun's motion are consistent with zero. 

any GD-1 member stars specifically. Therefore we have 
to search through the existing SEGUE spectra to iden- 
tify likely, or possible, member stars by position on the 
sky, CMD position and proper motion. In the previ- 
ous section we described that we used the ratio of the 
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Fig. 9. — Variation of the proper motions along the stream (cor- 
rected for the Solar reflex motion, assuming Vb = 220 km/s). The 
left panel shows the distribution of fi^ >1 for the stream candidate 
stars as a function of <f>i. The right panel shows the distribution of 
/j, c j >1 of the field stars selected using the same color-magnitude crite- 
ria as the stream stars. The variation of proper motions of stream 
stars with <pi is clearly visible in the left panel. Near <f>\ ~ 0° 
the proper motion of the stream is around —5 mas/yr, while at 
4>i ~ —50° it is around —8 mas/yr (corresponds to ~ 300 km/s at 
8 kpc). 

To construct the 6-D phase space distribution of the 
stream, the radial velocities are the remaining datum. 
By necessity the actual sample for which spectra, and 
hence radial velocities, will be available, will differ from 
the photometric sample just described. In this section 
we will use both the data from the SDSS/SEGUE survey 
(|Yannv et alj |2009) as well as radial velocities obtained 
by us with the TWIN spectrograph on Calar Alto, specif- 
ically targeting likely stream member stars. 

3.2.1. SDSS radial velocities 

SEGUE and SDSS only provide sparse spatial sam- 
pling of high latitude stars. SEGUE did not target 



stream/background probabilities — p BG ( M _ g!g _ r , r _^ r) 
select high probability members of the stream. Now 
additionally to that we also select the stars within the 
fj,a,Hs box (see Fig. [8]). That allows us to have a sample 
of stream stars with much less background contamina- 
tion, although the overall size of that sample is signifi- 
cantly smaller, since the the SDSS/USNO-B1.0 measure- 
ments of the proper mot ions were done for stars with 
r < 20 (|Munn et al.ll2004D . To illustrate how good the 
proper motion selection is when combined with the color- 
magnitude selection we show map of high probability 
stream member stars on Figure 1101 The stream is now 
clearly seen in individual stars. In Figure [10] we also 
overplot the location of existing SEGUE DR7 pointings, 
some of which cover the stream. Therefore we may ex- 
pect to find some stream members among the SEGUE 
targets in these fields. 

Figure [ED shows the SDSS/SEGUE radial velocity 
distribution as a function of <pi for those stars whose 
proper motions and color-magnitude position are con- 
sistent with stream membership, and which are located 
within 3 degrees from the center of the stream. The 
typical uncertainty of the SDSS/SEGUE radial veloci- 
ties is ~20 km/s. The filled red circles in this Figure 
show the subset of stars located within 0°.3 from the 
center of the stream and therefore represent the subset 
with very high membership probability, while the open 
black circles represent (spatially selected, \<j>2\ > 0.3°) 
background stars with similar proper motion and color 
magnitude. The filled symbols in Fig. [TT] clearly de- 
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lineate the radial velocity of the stream. Clumps of 
red circles are visible at ((pi, K a d)~(-25°,-100 km/s), (- 
30° ,-80 km/s), (-47°,0 km/s), (-55°,40 km/s). In or- 
der to perform the formal measurements of the radial 
velocities we performed a maximum likelihood fit by a 
model, consisting from two Gaussians one (wide) Gaus- 
sian was representing the background distribution, while 
the second (narrow) was modeling velocity distribution 
of stream stars. This fit gave us the following re- 
sults: (fa, V; od )=(-56 Q ,39±14 km/s), (-47°,-7±10 km/s), 
(-28°,-61±6 km/s), (-24°, -83±9 km/s). 

3.2.2. Color Alto radial velocities 

Since the SDSS/SEGUE radial velocities only provide 
constraints at a few points of the stream, we decided to 
obtain additional radial velocity information with tar- 
geted observations. We based the target selection for 
likely stream members on all the available information 
discussed in the previous sections: position on the sky, 
color-magnitude location and proper motions. We se- 
lected 34 stars likely members with r < 19 for the ob- 
servations. All these stars are within the sample plotted 
in the Figure 1101 These stars are mostly main-sequence 
and main-sequence turn-off stars with 18 < r < 19. 

The observations were performed the TWIN spectro- 
graph on the 3.5m telescope at Calar Alto observatory, 
during several nights of service observing in February 
2009. The TWIN spectrograph is the intermediate reso- 
lution spectrograph installed in the Cassegrain focus of 
the telescope. It consists from two separate spectroscopic 
channels (blue and red) behind the common entrance slit. 
The light from the slit is splitted into "red" and "blue" 
beams by a dichroic mirror. We used the blue and red 
arms of the spectrograph at a resolution of 4000 — 5000 
to observe the H^, Mgb lines and Call near-IR triplet 
respectively. The standard data reduction steps were 
applied to the dataset using custom written routines in 
Python language. The median signal to noise ratio per 
pixel was 7 for the blue spectra and 3 for red spectra. 

We used both the blue and the red spectra to com- 
pute the radial velocity of each star. The radial ve- 
locity of each star was derived by minimizing % 2 as a 
function of velocity shift of the template convolved with 
the appropriate Line Spread Function (LSF). The % 2 
for each star was a sum of the \ 2 f° r the blue and 
the red part of the spectra. As template in the blue 
spectral r ange we used the sp ectra from the ELODIE 
database (jPrugniel et al.l [20071 ) for stars of similar color 
and magnitude to the targeted ones and with low metal- 
licity [Fe/H] ~ — 2. In the red spectral range the tem- 
plate spectra was simply consisting from three lines of Ca 
triplet at 8498.02A, 8542.09A and 8662.14A. The error 
of each velocity measurement was determined using the 
condition A(x (V)) = 1- 

The Calar Alto measurements of the velocities to- 
gether with their errors are overplotted in blue sym- 
bols in Fig. [TTJ It is apparent from Fig. [TTJ that for 
—50° <</>i< —10° the targeting strategy was very suc- 
cessful, nicely delineating the projected velocity gradient 
along the stream. Overall out of 34 observed stars ~ 24 
stars belong to the stream and ~ 5 didn't have enough 
S/N for the velocity determination. Unfortunately the 
targeting near fa w +5° failed to identify stream mem- 



bers, probably because the stream there is less intense 
and further away. 

4. MODELING 

In the previous section we described the derivation of 
different stream properties such as distance, position on 
the sky, proper motion separately. In this Section we 
will map the stream in 6-D position-velocity space in a 
more co nsistent way, usin g all the available information 
(see e.g. iCole et alTl2008L for the application of similar, 
although simpler method to Sgr stream). This will pro- 
vide us with a set of orbit constraints along different 
sections of the GD-1 stream, which we will then model 
by an orbit to derive potential constraints. 

4.1. Positions on the sky and distances to the stream 

We start by characterizing the projected stream po- 
sition and its distance from the Sun through a maxi- 
mum likelihood estimate for a parametrized model of 
the stream P s tream(f, g — r,fa,fa) that describes it in 
4-dimensional space of photometric observables r,g — 
r,fa,fa: 



P s tream(r,g-r,fa,fa) = CMD(r, g - r, D(fa)) x 

1 ( (fc - fafi(fa)?\ 



Xl(fa) X 



2n(Tfa(fa) 



exp 



2<(0i) 



(2) 



Here fa,o(fa) is the fa position of the stream center 
on the sky as a function of (f>i, <70 2 (</>i) is the pro- 
jected width of the stream in fa, I (fa) is the "intensity" 
(i.e. the number density) of the stream as a function 
of fa, and D((f>i) is the distance to the stream. Fur- 
ther, CMD(r,g — r,D(fa)) is the normalized Hess dia- 
gram (i.e. the probability distribution in CMD space) ex- 
pected for the stream's stellar population at a distance of 
D(fa) after accounting for the observational errors. We 
construct that CMD based on the age and metallicity ob- 
tained in Section[2]and the isochrones from Girard fet al.l 



iMarigo et al.l (|20081 ) (assuming that 



and —Q)- In this model, P stream depends on four 

functions - I (fa), fa,o(fa), &fa(fa), D(fa) - which we 
take to be piecewise constant; i.e. for intervals 5 fa they 
simply become four parameters. 

For the field stars, Pbg( t , 9 — r , fa, fa) the analogous 
4D distribution is 

PBG{r,g~r,fa,fa) = I BG (fa,fa)xCMD B G(r,g-r,fa) 

, where lBc(fa,fa) is the 2D number density dis- 
tribution of the field stars around the stream and 
CMDbg( t , 9 ~ r i fa) the corresponding color- magnitude 
diagram. These functions are determined empirically 
from the data in parts of the sky adjacent to the stream 
(1021 ^ 0.5°). lBc(fa,fa) is determined by fitting the 
density of the stars in the fa, fa space by a polynomial. 
CM Dbg ( r , 9 ~ r , fa ) is determined by constructing the 
Hess diagrams using all the stars with 0.3° < \ fa\ < 5° 
in several fa bins. 

To simplify the determination of Pbg( t ,9 ~ r ,4>i,fa) 
and P s tream(r, 9~ r, fa, fa), we split the stream in several 
fa pieces, and consider I (fa), a^, D(fa) and fafi(fa) as 
constants within them. The log-likelihood for the mix- 
ture of the Pstream and Pbg distribution can be written 
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Fig. 10. — Stars that match the expectations for stream members with regards to proper motions (—15 mas/yr< fj,j, < —5 mas/yr, 
M</>2 I < 3 mas/yr), colors and magnitudes (but no rp2 filter), used in the candidate selection for radial velocity measurements. The stream 
can be clearly seen in distribution of individual stars. The locations of the SEGUE DR7 pointings are shown by gray circles. 

stars, as it varies noticeably along the stream (see Fig. [7]) 
and the reason of these variations is not clear. It is ap- 
parent from Fig. \TU\ that the projected stream position 
is very well defined, and that a distance gradient exists 
along the stream. 
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Fig. 11. — Radial velocities of likely stream stars (filled cir- 
cles). The Figure shows the radial velocities drawn from the 
SDSS/SEGUE (red circles) and Calar Alto spectra (blue circles). 
The red and blue symbols reflect the radial velocities of all stars 
matching in color, magnitude and proper motion that have posi- 
tions with 1 02 1 < 0.3°. The open circles represent SDSS/SEGUE 
velocities of similar stars but with \<j>%\ > 0.3°. The radial com- 
ponent of the Solar reflex motion (taking Vb = 220 km/s) was 
subtracted from all datapoints. 

as (here for convenience we introduce a as a fraction of 
stream stars instead of I(4>i)): 

ln(£) = 2 ln(a P s tream{n, 9i - n,<f>l,i,<f>2,i) + 

stars 

(1 - a)P B G(ri,gi - n,<pi til (f>2,t)) (3) 

and should be maximized with respect to the pa- 
rameters (ct0 2 , D(<pi), 02,o(0:l ) an d a). The max- 
imizatio n is perfor med using the Truncated Newton 
method (|Nash |[l984l) . The parameter errors are obtained 
using numerically com puted Fisher information (see e.g. 
ICox fc Hincklevlll^l . 

The left and central panels of FigurelT^lshow the result- 
ing estimates of the projected position and the distance 
of the stream, the parameters used in the subsequent or- 
bit fitting. We do not use the number density of stream 



4.2. Proper motions 

The likelihood maximization just described also results 
in stream membership probabilities for any given star i, 
Pstreumin, 9i - n,4>i,h 02,i)- This information can then 
be used to estimate via maximum likelihood the mean 
proper motions of different stream pieces, thereby ex- 
tending the observational estimates to the full 6-D space 
(r,g - r,0i,0 2 ,^ 1 ,/i0 2 ). 



Pstream,»(r, g - r, 0i , 2 , fifa , fifa ) 
P stream (l", 5 - r, 01 , 2 ) X y 

2iraf l 

(/% ~ M0i,o(0i)) 2 + 0<fe ~ M<fe.o(0i)) S 



exp 



(4) 



where we simply take the previously determined 
Pstream as a prior, to be modified by the Gaussian distri- 
bution in the 2D proper motions space. Here, the proper 
motion distribution is characterized by three functions, 
M0i,o(0i)> A*02,o(0i) and c^, which again we take to be 
piecewise constant. The distribution of the background 
stars in the 6-D space P B a^{r,g - r, 01,02,//^, (J,fa) is 
obtained empirically by binning the observational data. 
Then we construct again the logarithm of likelihood, con- 
sidering variations in four parameters: the number of 
stream stars, the proper motion of the stream in 4>i and 
02 and the proper motion spread ct m . This likelihood is 
then maximized and we determine the (ifa , [ifa and a M 
for different stream pieces. The right panel of the Fig- 
ure [T2] shows the resulting proper motion estimates as 
a function of 0i. These proper motions have not been 
corrected here for the Sun's reflex motion, which we will 
model in Section [5J 

Overall, the analysis presented in the previous sections 
has resulted in the best and most extensive set of 6-D 
phase-space coordinate map for a cold stream of stars in 
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TABLE 1 

Radial velocities from the Calar Alto Observations 





Star 


01 


02 


Vrad 






[deg] 


[deg] 


[km/s] 


SDSS 


TAfl A 1 f\K or 1 old 1 1 C 

JU941U5.oo+ololll.o 


—45.23 


—0.04 


o id o I c n 

lo.o ± 6.9 


bJJbb 


Tnri/i 7nc! ofi i ooonon o 
JU94;Uo.zo+ooz9o9.o 


—43.17 


—0.09 


on q _i_ i n i 


bJJbb 




—39.54 


—0.07 


2.9 dr. 8.7 


SDSS 


J095910. 43+363206. 6 


—39.25 


—0.22 


—5.2 dr. 6.5 


SDSS 


T 1 nnO OO n 1 i O ^7 A 1 -i o o 

J 100222. 01+3741 13. 3 


—37.95 


0.00 


1.1 dr 5.6 


SDSS 


Ti f\n<ion no i o ^7 i n , i n o 

J1UUzz2.(Jz-t-o74U49.z 


—37.96 


—0.00 


— 11.7 dr. 11.2 


SDSS 


ti m noo no i onoonn o 

J 101033.02+393300.8 


—35.49 


—0.05 


—50.4 dr 5.2 


SDSS 


Timiinno i on /i /i tro n 

J 101 110. 08+394453. 9 


—35.27 


—0.02 


on n I i o o 

—30.9 ± 12.8 


SDSS 


Timor/ oo i onccnir 

J 101254. 83+395525.6 


—34.92 


—0.15 


—35.3 ± 7.5 


SDSS 


T1 n 1 <ii o n r i i aat 1 1 o o 

J 101312.05+400613.3 


—34.74 


—0.08 


on n 1 n o 

—30.9 ± 9.2 


SDSS 


Ti n 1 ^no i r i a r\ a ^7 a o 

J 101702. 15+404747. 3 


—33.74 


—0.18 


r-i a o i no 

—74.3 ± 9.8 


SDSS 


J101951. 76+412701. 5 


-32.90 


-0.15 


-71.5 + 9.6 


SDSS 


J102216.20+415534.7 


-32.25 


-0.17 


-71.5 + 9.2 


SDSS 


J103003.87+434351.7 


-29.95 


-0.00 


-92.7 + 8.7 


SDSS 


J104341.92+460224.7 


-26.61 


-0.11 


-114.2 ± 7.3 


SDSS 


J104840.98+464922.1 


-25.45 


-0.14 


-67.8 + 7.1 


SDSS 


J105036.96+472000.1 


-24.86 


0.01 


-111.2 ± 17.8 


SDSS 


J110711. 27+494415. 9 


-21.21 


-0.02 


-144.4 + 10.5 


SDSS 


J114242.08+533841.4 


-14.47 


-0.15 


-179.0 ± 10.0 


SDSS 


J114724.59+535546.8 


-13.73 


-0.28 


-191.4 + 7.5 


SDSS 


J115116.08+542142.7 


-13.02 


-0.21 


-162.9 + 9.6 


SDSS 


J115326.06+542930.6 


-12.68 


-0.26 


-217.2 ± 10.7 


SDSS 


J115404.06+543511.4 


-12.55 


-0.23 


-172.2 + 6.6 



TABLE 2 

Stream positions 



01 


02 


[deg] 


[deg] 


-60.00 


-0.64 + 0.15 


-56.00 


-0.89 + 0.27 


-54.00 


-0.45 + 0.15 


-48.00 


-0.08 + 0.13 


-44.00 


0.01+0.14 


-40.00 


-0.00 + 0.09 


-36.00 


0.04 + 0.10 


-34.00 


0.06 + 0.13 


-32.00 


0.04 + 0.06 


-30.00 


0.08 + 0.10 


-28.00 


0.03 + 0.12 


-24.00 


0.06 + 0.05 


-22.00 


0.06 + 0.13 


-18.00 


-0.05 + 0.11 


-12.00 


-0.29 + 0.16 


-2.00 


-0.87 + 0.07 



TABLE 3 

Stream distances 



01 


Distance 


[deg] 


[kpc] 


-55.00 


7.20 + 0.30 


-45.00 


7.59 ± 0.40 


-35.00 


7.83 ± 0.30 


-25.00 


8.69 + 0.40 


-15.00 


8.91 ± 0.40 


0.00 


9.86 + 0.50 



the Milky Way. 

5. ORBIT FITTING 

If we can assume that all the stream stars lie close to 
one single test-particle orbit, then our phase space map 
of the GD-1 stream should not only define this orbit, but 
at the same time constrain the Milky Way's potential. In 



TABLE 4 

Stream proper motions 



01 
[deg] 


[mas /yr] 


^02 

[mas /yr] 


[mas / yr] 


-55.00 


-13.60 


-5.70 


1.30 


-45.00 


-13.10 


-3.30 


0.70 


-35.00 


-12.20 


-3.10 


1.00 


-25.00 


-12.60 


-2.70 


1.40 


-15.00 


-10.80 


-2.80 


1.00 



general the stars in the tidal stream have slightly differ- 
ent energies and angular momenta, but the assumption 
that the stream stars are moving along the same orbit is 
plausib le , especially if the s tream is thin and is near peri- 
center (jDehnen et al.l [20041. & private communication). 
But it is not straightforward to quantify the quality of 
such an approximation. For now we simply fit an or- 
bit to our 6-D map of available observational data: the 
position on the sky, 4>2{4>i)i proper motion /2(</>i), dis- 
tance to the stream D(<fii) and radial velocity V ra d(4>i)- 
For each assumed potential, we will determine the best 
fit orbit, but then marginalize over the orbits to deter- 
mine the range of viable gravitat ional potentials. This 
analysi s ext ends earlier ef f orts b y iGrillmair fc Dionatosl 
(|2006bf ) and lWillett et al.l (|2009f ) who have presented or- 
bit solutions for GD-1 . However, we can now draw on 
a much more extensive set of observational constraints. 
We also explore the fit degeneracies. Given that our 6-D 
phase-space map of the GD-1 stream spans only a lim- 
ited range in R and z (as seen from the Galactic center 
11 kpc<R< 14 kpc, 5 kpc<z< 9 kpc) it proved useful 
to consider very simple parametrized potentials at first. 
Further it proved necessary to consider what prior infor- 
mation we have on the Sun's (i.e. the observers) position 
and motion, as well as on our Milky Way's stellar disk 
mass. 

5.1. One component potential 

The stream is located at Galactocentric (R, z) « 
(12,6) kpc, a regime where presumably both the stel- 
lar disk and the dark halo contribute to the potential , 
and its flattening. Of course, the stream dynamics are 
solely determined by the total potential, and therefore 
we consider first a simple single-component potential, the 
flattened logarithmic potential 

*(x,y,z)='£ln(x 2 +y 2 + (j^jj , (5) 

which has only two parameters: the circular velocity V c 
and the flattening g$. Note that (1 — qderisitv) ~ 3(1 — a®) 
for m oderate flattening (e.g. p. 48 of IBinnev fc Tremaind 
1987). Such a simple potential seems justified as the 
stream stars are only probing a relatively small range in 
R and z. 

In practice, we fit an orbit to the 6-D stream map 
by considering a set of trial starting points in the 

Galaxy, specified by the initial conditions (X(0),X(0)) 
in standard Cartesian Galactic coordinates. Together 
with an assumed gravitational potential this predicts 

(X(t),X(t)), which can be converted to the observables, 
4>2(4>i), /3(0i), D{<f>i), V ra d(4>i) and then compared to the 
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Fig. 12. — Summary of photometrically derived stream properties based on maximum likelihood fits to chunks of the stream, drawing 
on SDSS photometry and astrometry (see Section 0}. The left panel shows the positions of the stream on the sky in <f>i,<p2 coordinates. 
The middle panel show the measurements of the distances as a function of <j>\. The right panel shows the statistical proper motions of the 
stream stars (without the correction for the Solar motion); red circles show the )i r f >l (the proper motion along the stream), blue squares 
show the 11^2 (the proper motion perpendicular to the stream). 



6-D observations (Fig. [13]). For each [(X(0), X(0)\$(X)] 
we can evaluate the quality of the fit by calculating % 2 , 
summing over all data points, shown in Fig. 1131 For 
any given $>(X), \ 2 can be then minimized with respect 

to the orbit, i.e. (X (0) , X (0)) , providing the 'best fit' 
orbit in this potential and the plausibility of that poten- 
tial. The minimization is performed using the MPFIT 
code ( M arkwardtJ |2009|) implementing the Levenberg- 
Marquardt technique (jMarquardtl 119631 ) translated into 
Python 7 . The data used to constrain the potential are 
given in the Tables [TJ [3J [31 H] (except the SDSS measure- 
ments of the radial velocities which are given in the end 
of Section [3]) . 

It is crucial to note that the conversion of (X(t) 7 X(t)) 
to the space of observables depends on the position and 
motion of the observer, i.e. on distance from the Sun 
to the Galactic center (i? ) and on the 3-D velocity of 

the Sun in the Galaxy rest- frame (Vq). At this stage we 
adopt i?n = 8-5 kpc based on recent de terminations (e.g. 
iGhez et al.ll2008b iGillessen et al.ll2009f ). but later we will 

relax this. The second parameter Vq = Vlsr + AVlsr 
(where Vlsr. is the velocity of the Local Standard of 
Rest and AVlsr is the Sun's velocity relative to the 
LSR) is linked to the fitting not only through conver- 
sion of the observable relative stream velocities to the 
GC rest system, but also conceptually through the plau- 
sible demand that &(X) and in particular V c (Ro,0) also 
reproduces Vlsr- In this way, constraints on the poten- 
tial flattening can be derived by considering in the 

disk plane (Vlsr) and the plane of the GD-1 stream. 
The velocity of the Sun relative to the Local Standard 

of Rest (LSR) AVlsr is qu ite well known from th e 
HIPPARCOS measurements (jDehnen fc Binnevl [l998h : 

AV L SR[km/s] = I0e x + 5.25e y + 7.l7e z . The velocity 
of the LSR, i.e. V c (Ro, 0) has a considerable uncertainty 

7 http: //code . google . com/astrolibpy/ 



(IBrand fc Blitz! [1991 IGhez et all [20081: IXue et all [20081: 
iReid et al.ll2009t ). Initially we will consider the velocity of 
the LSR simply a consequence of the assumed potential, 
i.e. V LS R = V c {R o ,0). 

Figure [T3] illustrates the result of such fitting, by over- 
plotting the best fit orbit for the plausible potential with 
V c = 220 km/s and g$ = 0.9 over observational data. 
It is clear that even for the simple flattened logarithmic 
potential, an orbit can be found that reproduces the ob- 
servables well. This fit and Figure serve to illustrate a 
few generic points that also hold for orbit fits in differ- 
ing potentials: the stream moves on a retrograde orbit 
and it is near pericenter, wh ere the stream is ex pected 
to approximate an orbit well (jDehnen et a l. 2004). After 
fitting a first orbit, we may also note its global parame- 
ters (see Fig. [T71 for a 3-D map of the orbit): pericenter 
is at 14 kpc from the GC; apocenter is at 26 kpc; and 
the inclination is 39°. 

For any given potential $(X\V c ,q®) we can find best- fit 

orbit by marginalizing over (X(0),X(0)) to see what our 
6-D map of GD-1 implies about the relative plausibility 
of different V c and g$: Figure [T4l shows the log-likelihood 
surface for the potential parameters (V c ,q$)] note again 
that this fit neglects all other prior information on V c 
at the Sun's position. The contours show ler, 2er, 3c 
confidence regions on the parameters, derived from the 
S(\n(£)) values for tw o degrees of freedom (i.e. a two 
parameter fit)(see e.g. lLampton et aT]|1976D . The insets 
at the left and bottom show the marginalized distribu- 
tions for single parameters. The best fit values with the 
2-sided 68% confidence intervals are V c = 221±H km/s 

and = 0.87t{j;^. Fi gure [T4lillustrates that the flatten- 
ing parameter g$ is quite covariant with the equatorial 
circular velocity V c . An extreme example may serve to 
explain this covariance qualitatively. If the stream went 
right over the pole (z-axis), then the local force gradient 
would be proportional to V c x g$ (Eq. [5]) . Information 
about the potential flattening must therefore come from 
combining kinematics and dynamics in the disk plane 
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Fig. 13. — The data-model comparison for the best fit orbit in a flattened logarithmic potential (Eg. 151 with V c = 220 km/s and g$=0.9. 
The color data points with error bars shows the observational data, while the black lines show the model predictions for the orbit with 

X(0) = (-3.41, 13.00, 9.58) kpc, X(0) = (-200.4, -162.6, 13.9) km/s. The top left panel shows the positions on the sky, the top right panel 
shows the proper motions, the bottom left panel shows the distances, the bottom right panel shows the radial velocities. On the top right 
panel, red circles and thin line show /i^ , while blue squares and thick line show fi^ . 



with the information from GD-1 . 

The fit shown in Fig. [14] asks the data not only to con- 
strain the potential at the stream location and determine 
the stream orbit, but also to infer the Sun's motion (or 
at least Vlsr) from its reflex effect on the data. Clearly 
providing a prior on V c (Rq, 0) is sensible, especially if we 
care about constraints on the potential flattening. We 
consider the constraints that arise from the Sun's reflex 
motion with respect to the Galactic c enter the most ro - 
bust and sensible prior in this context. IGhez et al.l ((2008) 



recently combined radio data (jReid fc Brunthalerl 120041) 
with near-IR data on the Galactic center kinematics to 
arrive at V c (Rq,0) = 229 ± 18 km/s. It is also notice- 
able that our own constraint on V c (Ro,0) from Fig. Q3] 
221I20 km/s, is close both in value and uncertainty to 
the estimate of IGhez et al.l (|2008l ). which is based on a 
completely disjoint dataset and approach. 

Fig. [T5l shows the resulting log-likelihood contours and 
la, 2a, 3<7 confidence regions after applying this prior on 
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Fig. 14. — The log-likclihood surface of the orbit fit for the family 
of flattened logarithmic potentials (Eq. with different circular 
velocities V c and Battenings g$ with a flat prior on V c . Note that 
V c enters both into the model velocities of the stream stars and into 
the correction of all three velocity components for the Sun's motion. 
The contours show the la, 2a and 3c confidence regions. The inset 
panels at the bottom and on the left show the ID marginalized 
posterior probability distributions for V c , q& respectively. The gray 
line in t he bottom panel s hows the probability distribution for the 
V c from lGhez eTa l. ( 2008), which we shall use as a prior in Fig. 1151 



log(£) 




160 180 200 220 240 260 280 300 
V c . [km/s] 



Fig. 15. — The log-likclihood surface of the orbit fit for the fam- 
ily of flattened logarithmic potentials (Eq.[5]l with different circular 
velocities V c and flattcnin gs q§, but no w with a prior on the V c 
of 229±18 km/s from Gh ez et al.l II2008T ). The likelihood was also 
margin alize d over the Gaussian prior on Ro = 8.4 ± 0.4 kpc. As 
on Fig. 1141 the contours show the la, 2a and 3a confidence re- 
gions. The inset panels at the bottom and on the left show the ID 
marginalized posterior probability distributions for V c , q& respec- 
tively. The gray line in the bottom pan e l show s the adopted prior 
distribution for the V c from lGhez et al.l 12008T ) 

the V c (229±18 km/s). Note that likelihood on Fig.[T5lis 
also marginal ized over Rg with a Gaussian prior (R = 
8.4 ±0.4 koc lGhez et Zll2008h . 

Fig. [15] shows that the posterior probability distribu- 
tion on V c has slightly changed to V c (Rn) = 224^^ km/s 
with noticeably smaller error bar compared to the value 
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Fig. 16. — The data-model comparison for a set of best-fit or- 
bits in different logarithmic potentials (Eq. [5J with three differ- 
ent (V c ,g$) parameters values (180 km/s, 1.1), (220 km/s, 0.9), 
(260 km/s, 0.8.) The colored data points with error bars show the 
observational data, while the black lines show the model predic- 
tions (different line styles show the orbit models in different poten- 
tials). The top left panel shows the positions on the sky, the top 
right panel shows the proper motions, the bottom left panel shows 
the distances, the bottom right panel shows the radial velocities. 
On the top right panel, red circles and thin lines show , while 
blue squares and thick lines show fi,p 2 . 

from iGhez et al.l (|2008l) . Fig US] also shows us the slight 
improvement comparing to Fig. [14] of flattening con- 
straints: g<j, = 0.87^q This means that the total po- 
tential appears to be oblate (in the radial range probed); 
this may not be surprising, as the stellar disk - which is 
manifestly very flattened - contributes to the total po- 
tential. 

Fig. [161 illustrates how well the best-fit orbits for differ- 
ent potentials (&(X\V c ,q$) can mimic one another other 
in the space of observables. This is the source of the 
parameter covariances shown in Fig. [14] and 1151 

The fitting of the orbit shown in Fig. [T3] allows us to 
make an estimate of the line-of-sight velocity dispersion 
in the stream, by comparing the dispersion of the ra- 
dial velocity residuals with the accuracy of individual 
radial velocities. Assuming that without observational 
uncertainties the residuals from the orbit fit should be 
Gaussian distributed with zero mean we performed the 
maximum likelihood fitting of the residuals, which gives 
a 90% confidence upper limit on the velocity dispersion 
of stars in the stream of ~ 3 km/s. 

Before we aim at separating possible flattening contri- 
butions from the halo and disk, it is worth commenting 
on the accuracy and limitations of our estimate of g$. 
In the range (R, z) ~ (12, 6) kpc no other direct obser- 
vational constraints on the potential shape exist in the 
literature, and hence our estimate of g$ = 0.87^q'q4 is a 
new and important contribution. On the other hand, an 
error of 5q$ ~ 0.05, especially when translated into the 
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Fig. 17. — 2-D projections of the orbit in the Galactic rectangular 
coordinates. The position of the Sun is shown by a red circle. The 
Galaxy is shown by a cloud of points and the gray arrow shows the 
direction of the galactic rotation. The black arrow shows the direc- 
tion of the orbital movement of the stream stars. The orbit is the 
best fit orbit for the V c = 220 km/s, q$ = 0.9 logarithmic potcn- 
tial(Eq. [5}. The orbit for the 3-component potentiaKEqns, l6l 171 l8t 
is almost undistinguishable from the orbit in logarithmic potential. 

flattening error of the equivalent scale-free mass distri- 
bution, may not appear as particularly helpful in model 
discrimination, or as impressively accurate. Especially 
as a manifestly cold stellar stream spanning over 60° on 
the sky may seem ideal for mapping the potential at first 
glance. 

5.2. Constraints on the shape of the dark matter halo 
from a bulge, disk, halo 3-component potential 

In the previous section we constrained the parame- 
ters of a simplified MW potential, the spheroidal log- 
arithmic potential. It is clear that the MW potential 
at the position of the stream must depend explicitly 
on the sum of baryonic Galaxy components (bulge and 
disk) and on the dark matter halo. We now explore 
whether our constraint on the shape of the overall po- 
tential, q$ ~ 0.9, permits interesting statements about 
the shape of the DM potential itself. At the distance 
of (R, z) ps (12, 6) kpc the contribution of the disk to 
the potential should still be relatively large, weakening 
or at least complicating inferences on the shape of the 
DM distribution. 

We adopt a three-component model of the Galaxy po- 



tential, choosing one tha t is widely used in the mod- 
eling of the Sgr stre am (|Helmil 120041: lLaw et al.l 120051: 
iFellhauer et al.l |2006T > and reproduces the galactic rota- 
tion curve reasonably well. 

The model consists of a halo, represented by the loga- 
rithmic potential The numeri cal values of par a meter s of 
the potential were taken from IFellhauer et all (|2006t ) 

<t>haio(x,y,z) = %Mn ( x 2 + y 2 + ( —?—) + d 2 j , 

* \ \q<S>,halo/ J 

.(6) 

where we have adopted d = 12 kpc from the previous 
authors. The disk is represented b y a Miyamoto-Nagai 
potential (jMivamoto &: Nagailll975[ ). 

^disk(x, y, z) = = (7) 

Jx 2 +y 2 + (b + \Jz 2 + c 2 ) 2 

with b — 6.5 kpc, c = 0.26 kpc. The bulge is modeled as 
a Hernquist potential: 

. . GMb 

®bulge(X,y,Z) = ■ (8) 

r + a 

with M b = 3.4 x 10 10 Af Q , a = 0.7 kpc 

As in the previous section, for any given set of po- 
tential parameters, we can find the best-fitting stream 
orbit and compute y 2 of the fit. In the current paper 
we do not make any attempts to fully fit all the pa- 
rameters of the MW potential, but we try to make an 
estimate of the MW DM halo flattening. We take the 3- 
componcnt potential and fix all but 3 parameters — disk 
mass Mdisk, circular velocity of the halo Vhaio and the 
flattening of the halo <7$,/ la /o- On a 3-D grid of these pa- 
rameters we perform a % 2 fit. Figure [TBI shows the results 
of such fit after marginalization over the orbital parame- 
ters (A(0), A(0)), circ ular velocity of th e halo v^aio with 
a Gaussian prior from lXue et al.l (|2008T > and a Gaussian 
prior on the circu lar velocity at the Sun's radius from 
iGhez et all (|2008f ). The Figure clearly illustrates that in 
the case of the 3-component potential, the current data 
is unable to give a significant new insights on the flat- 
tening of the DM halo. We can only say that at 90% 
confidence q<s>,haio > 0.89. We note however that for the 
future analysis, if a multi-component model for the po- 
tential is used then more prior information is required, 
i.e. not only on Vlsb. & n d Vhaio but on M^sk and other 
parameters of the potential. 

6. DISCUSSION AND CONCLUSIONS 

In this paper we have presented a thorough analysis 
the GD-1 stream combining the publicly available SDSS 
and SEGUE data with follow-up spectroscopic observa- 
tions from Calar Alto. The combination of the photomet- 
ric SDSS observations, USNO-B/SDSS proper motions, 
SDSS, SEGUE and Calar Alto radial velocities allowed 
us to construct a unprecedented 6-D phase-space map of 
this very faint (29 mag/sq. arcsec) tidal stream. The 6-D 
phase-space map of the stream, spanning more than 60° 
on the sky, provided the opportunity not only to fit the 
orbrt as lGrillmair fc Dionatosl (|2006b[ ) and lWillett et al.l 
(2009) have done previously but to explore what con- 
straints can be placed on the MW potential. 
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Fig. 18. — The log-likelihood surface of the orbit fit for a 3- 
component potential (Eq. [6] [7] and [SJ with different disk masses 
Mdiski halo circular velocities v^alo and halo Battenings q$ y halo- 
The likelihood was marginalized over the circular velocity of the 
halo Vhalo with the Gaussian prior on v^ a i = 170 ± 15 km/s from 
IXue et al. ( 2008) and th e Gaussian prior on Vlsr = 229±18 km/s) 
from[Ghcz ct al. (2008) . The contours show the lcr, 2<r and 3<r 
confidence regions. The inset panels at the bottom and on the left 
show the ID marginalized posterior probability distributions for 
M disk , q<s,,halo respectively. 

The analysis is based on the assumption that the 
stream stars occupy one orbit. In detail, of course, dif- 
ferent stars on the stream have slightly different values 
of conserved quantities and therefore lie on slightly dif- 
ferent orbits. Effectively, our analysis depends on these 
being small when compared with an orbital uncertain- 
ties. The magnitude of the departure of the stream from 
a single orbit will, in detail, be a function of the projcn- 
itor and the distruption process; as these details became 
understood the model can be refined. 

We found that the distribution of stream stars in phase 
space can be well fit by an inclined eccentric orbit in the 
spheroidal logarithmic potential. After marginalization 
over the stream orbital parameters we derive a circular 
velocity V c = 221±2o km/s and flattening = 0.871^3. 
This measurement has been made without the use of any 
information other than that in the GD-1 stream itself. 
It is important that the information available in the ob- 
servations of the stream is very sensitive to the V c , the 
circular velocity at the Sun's position. The reason for 
that is that the stream extends more then 60° on the sky 
and therefore both the radial velocities and the proper 
motions have components coming from the projection of 
the Sun's motion. 

If we combine our circular vel ocity measu r ement with 
existing prior on the V c from Gh ez et alJ (|2008| ) and 
also marginalize over the dis tance from the S un to 
the Galactic center using the IGhez et alJ ((2008) prior 
(Ro = 8.4 ± 0.4 kpc), we further tighten the error bar on 
V c — 2241 {J km/s and on the flattening of the potential 

q$, = 0.87l ; 04 . Our measurement of the V c is the best 
constraint to date on the circular velocity at the Sun's 
position, and the measurement of g$ is the only strong 
constraint on g$ at galactocentric radii near R ~ 15 kpc. 
The measurement of the flattening of the potential 



qq, = 0.87_o'o4 describes only the flattening of the over- 
all Galaxy potential at the stream's position (R, z) w 
(12, 6) kpc where the disk contribution to the poten- 
tial is presumably large. Unfortu nately the data o n 
the GD-1 s tream combined with the IGhez et al.l (|2008f ). 
IXue et all (|2008l ) priors on V c and Vhalo are not enough 
for separating the flattening of the halo from the flat- 
tening of the total Galaxy potential. So we are un- 
able to place strong constraints on flattening of the 
MW DM halo; we put a 90% confidence lower limit at 
q<s>,haio > 0.89. 

Despite the negative result on the measurement of the 
MW DM halo flattening, we note that the data from the 
GD-1 stream is able to give strong constraints on two 
important Galaxy parameters. We claim that that our 
dataset on the GD-1 stream should now be combined 
with other available MW kinematical data (from other 
stellar streams, BHBs, MW rotation curve) in order to 
tighten the existing constraints on the Galaxy param- 
eters. It is important that the constraints on Galaxy 
potential based on the GD-1 stream dataset are, to large 
extent, model- independent and purely kinematic i.e. the 
constraints on the V c come to large extent from the pro- 
jection effects and manifest themselves in proper motions 
and radial velocities. 

We think that it is quite surprising that such a long 
(60°) stream with full 6-D map didn't allow us to con- 
train large number of parameters of the MW DM halo 
and other Galactic components. We think that there 
are several reasons for that, but the most important is 
that while the observed part of the orbit spans ~ 70° 
on the sky from the Sun's point of view, the orbital 
phase spanned by the stars as seen from the Galactic 
center, is only ~ 40°. Furthermore, the observable part 
of the stream occupies the perigalacticon part of the or- 
bit, so the range of Galactocentric distances probed by 
the stream is small. That makes it plausible why orbits of 
different eccentricities, and hence of different azimuthal 
velocities at their pericenter, can match so closely the 
same set of 6-D coordinates. Since very cold streams take 
many orbits to spread a substa ntial fraction of 27r in or- 
bital phase (see e.g. also Pal 5: IOdenkirchen et alJl200H 
iGrillmair fc Dionatosll2006"al) . all future analyses of yet- 
to-be discovered streams will particularly need to con- 
sider the trade-off between the conceptual and practical 
attractiveness of 'cold' streams and the near-inevitable 
limitations of their phase coverage. 

In addition to the weakness coming from phase cov- 
erage, our analysis at this point must rely on photo- 
metric distance estimates; these have random errors of 
~1Q%, after an empirical distance correction to the best 
fit isochrones that is of the same magnitude (see Sec- 
tion^. While the proper motions that we derived for 
ensembles of stream stars are unprecedented for a stel- 
lar stream in the Milky Way's outer halo; yet, the corre- 
sponding velocity precision, especially when compounded 
by distance errors, is still the largest single source of un- 
certainty in the fitting (e.g. our tests have shown that 
overestimating heliocentric distances to the stream leads 
to the overestimated measurement of V c ). Both the defi- 
ciencies of the distances and proper motions on GD-1 will 
be largely alleviated after the launch of the Gaia satellite, 
which will allow much tighter constraints on the Galactic 
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potential. 

But before the launch of Gaia, we suggest that further 
observations of radial velocities of stars in the stream 
(with deeper spectroscopy) and improveme nts in proper 
motio n precision (with e.g. Pan-STARRS Kaise r et al.l 
2002) should be able to reduce the error bars on Galaxy 
parameters significantly. It is also important to calibrate 
properly the distance to the stream, which may be done 
by confirming several probable BHB candidates in the 
stream. 

We sugge st that any attempt t o fit the Galactic poten- 
tial(such as iWidrow et al.l 12008) now should not ignore 
the dataset on GD-1 and should incorporate it into their 
fits. 

The orbital parameters, which we have measured for 
the GD-1 stream, are mo re or less consistent with those 
from lWillett et all (|2009( ): pericenter is at 14 kpc from 
the GC, the apocenter is at 26 kpc, the orbit inclination 
with respect to the Galactic plane is 39°. We have 
also estimated the total stellar mass associated with the 
stream to be ~ 2 x 1O 4 M , which together with the rel- 
atively small stream width of <~ 20 pc suggests that the 
progenitor of the stream was a globular cluster, although 
we cannot completely rule out the dwarf galaxy progeni- 
tor. We have also determined the 90% confidence upper 
limit for the velocity dispersion of the stars in the stream 
~ 3 km/s which do not contradict significantly neither 
globular cluster nor dwarf galaxy progenitor hypothesis. 
Given the length of the observable part of the stream 
of ~ 10 kpc, and the velocity dispersion < 3 km/s the 
age of the stream can be estimated to be greater than 
1.5 Gyrs (assuming that the progenitor is in the mid- 
dle of the observed part of the stream). It is interesting 
that the stream managed to evade possible destruction 
by interaction w ith DM subhalos orbiting around MW 
(|Carlbersl I2009T) . Although the dumpiness observed in 
the stream may be attributed to these past interactions 
(Koposov et al 2010 in prep.). 

Overall in this paper we have illustrated a method of 
analyzing the thin stellar stream using all the available 
information on it and further utilizing that to constrain 
the Galaxy potential. We believe that in the epoch of 
Pan-STARRS, LSST and especially Gaia, which will give 



us a wealth of new information on the MW halo, stream- 
fitting like that presented in this paper will be extremely 
useful and productive. 
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APPENDIX 

TRANSFORMATION OF EQUATORIAL COORDINATES TO STREAM COORDINATES <j> u <f> 2 



cos(<K)cos((/>2)\ /-0.4776303088 -0.1738432154 0.8611897727\ /cos(a) cos(<5)\ 
sin(0i)cos(0 2 ) = 0.510844589 -0.8524449229 0.111245042 x I sin(a) cos(5) 
sin(0 2 ) / V 0.7147776536 0.4930681392 0.4959603976/ \ sin(tf) / 
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